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RESUMEN 

El resumen sera traducido al espanol por los editores. 
This article describes the production of primordial 4 He nuclei in an inhomogeneous 
universe. The baryon distribution is spherically symmetric and consists of a high 
density inner region and a low density outer region. As the temperature decreases 
neutrons diffuse to the outer region until they are homogeneously distributed, and 
protons may be redistributed depending on how fast diffusion occurs. Nucleosyn- 
thesis occurs earlier in the inner region and neutrons diffuse back to that region. 
The rapidity of diffusion determines how much 4 He is ultimately produced. 

ABSTRACT 

This article describes the production of primordial 4 He nuclei in an inhomo- 
geneous universe. The baryon distribution is spherically symmetric and consists of 
a high density inner region and a low density outer region. As the temperature 
decreases neutrons diffuse to the outer region until they are homogeneously dis- 
tributed, and protons may be redistributed depending on how fast diffusion occurs. 
Nucleosynthesis occurs earlier in the inner region and neutrons diffuse back to that 
region. The rapidity of diffusion determines how much 4 He is ultimately produced. 
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f. INTRODUCTION 

In 1999 the author wrote a big bang nucleosynthesis code corresponding to a universe with an inhomoge- 
neous baryon distribution. The code envisions the universe as a lattice of spheres and treats one sphere as 
characteristic of the rest of the universe. The sphere has an initial radius ( the distance scale ) of 25000 
cm. at temperature T — 100 GK. The sphere consists of a high density inner region and a low density outer 
region. The boundary radius r& between the regions is set to r$/2 and the density contrast ratio R p = 
800:1. The baryon-to-photon ratio rj in this model is 3 xl0~ 10 ( rjxo = 3 ). The sphere is divided into a core 
and 63 spherical shells. 

The code solve for the number density n(i, s) of isotope species i contained in spherical shell s. The evolution 
of n(i, s) obeys the equation (Kainulainen et al. 1998) 
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The first term accounts for nuclear and weak reactions occuring between isotopes within shell s. The second 
term corresponds to the expansion of the universe. And the last term accounts for diffusion of isotope i between 
shells, though in the model only neutrons diffuse. 
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Fig. 1. Number densities of neutrons, protons and 4x 4 He 
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Fig. 2. Neutron Reaction Rates ( shell 3 and shell 62 ) 
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2. RESULTS FOR Ri = 25000 CM, mo = 3 

Figure (1) shows the number densities of free neutrons ( n(n, s) ), protons ( n(p, s) ) and 4 x 4 He ( 4n( 4 He, s) 
). The horizontal axes correspond to the temperature T from 100 GK down to 0.01 GK and the distance scale 
ri normalized to 64 dimcnsionless units, and the vertical axis is number density n(i, s) normalized to the initial 
baryon density ( ribo ). At T = 100 GK number densities in the inner shells vs. the outer shells are large. 
Weak reactions convert neutrons to protons until they no longer become significant at around T = 13 GK. 
Figure (2a) shows the neutron-proton conversion rates for shell 3, a high density inner shell ( normalized to 
ribo and to the expansion rate of the universe a ). Figure (2a) also shows that around T = 6 GK neutrons 
start to diffuse significantly into the lower density outer shells. But since the conversion rates are no longer 
effective at this temperature the protons in the inner shells remain there. At around T — 2.5 GK the neutrons 
are homogeneously distributed throughout the whole sphere. 

Nucleosynthesis occurs at about T — 0.9 GK, when the nuclear reaction n + p <-> d + 7 falls out of 
equilibrium. Other nuclear reactions then build heavier nuclei, but most baryons wind up as part of He nuclei. 
Because of higher density nucleosynthesis occurs earlier in the inner shells, depleting neutrons. So neutrons 
diffuse back into the inner regions as shown on the right of Figure (2a) . Figure (2b) shows the same reaction 
rates for shell 62, a low density outer shell. The rate for n + p^d + 7is lower and flatter than in Figure 
(2a). Figure (1) shows that most neutrons back diffuse into the inner shells and then undergo nucleosynthesis, 
leading to 4 He production overwhelmingly concentrated in the inner shells. 

The final mass fraction X4 He of the whole model is 0.243664. This is the result for rj = 25000 cm and 
7710 = 3 Table (1) shown below lists Xi He for n = 100 cm up to 6.3 xlO 6 cm. For r, less than 25000 cm 
neutron diffusion occurs early enough to coincide at least partly with the time when the conversion reactions 
are significant. In the inner regions the reactions convert protons into neutrons that then diffuse to the outer 
regions where they get converted back to protons. But as gets larger more protons remain in the inner 
regions and A4 He increases. r» = 25000 cm is the turning point when all the protons in the innermost regions 
remain there. For larger the back diffusing neutrons can't reach the innermost regions before nucleosynthesis 
is completed. So with fewer neutrons in the inner regions X<t He decreases. At ri = 7.9 x 10 5 cm though diffusion 
becomes significant at T = 0.9 GK, right during nucleosynthesis. For larger r, nucleosynthesis occurs before 
neutron diffusion and the inner regions, with high proton and neutron densities, produce enough He to raise 
A4 He significantly. 



ri in cm. 


A4 He 


10 2 


0.239029657 


1.58489 xlO 3 


0.239369667 


2.51189 xlO 4 


0.243663794 


10 5 


0.241847338 


7.94328 xlO 5 


0.229781933 


6.30957 xlO 6 


0.250711334 



Table (1): 4 He, where i) = 3.0x lO" 10 . 

3. CONCLUSIONS AND FUTURE RESEARCH 

Figure (3a) shows A4 He for a range of r\ from 10~ 10 to 3 x 10~ 9 as well as the range of r, listed in Table 
(1) For all values of r\ shown one can see the turning points at which Xi He increases and then decreases and 
then increases again. The code can generate contour maps for 67 other isotope species as well. Comparing 
the contour maps of A4 He , Y(d)/Y(p) ( Y(d) is the abundance of deuterium ), and Y( 7 Li + 7 Be)/Y(p) to the 
2a ranges of recent observations ( IT, 1998; (Izotov & Thuan 1998) BT, 1998; (Buries k Tytlcr 1998) PWSN, 
1998 (Pinsonncault ct at. 1998) ). 



= 0.244 ±0.002 
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Fig. 3. A4 He and 10 4 x AX4 He 

F(d)/y(d) = (3.4 ±0.3) x 1CT 5 
log[y( 7 Li)/y(d)] = -9.45 ±0.20 

the author determined an observational range for rj from 4.7 xlO -10 to 6 xlO -10 . 

A planned article will describe 4 He production in greater detail, describing results for values of rj other 
than 25000 cm. R p , and the geometry of the model can be varied in the code. In the future the author 
will look at results with different values for these parameters. Also, the code has written into it subroutines to 
calculate the neutrino heating effect. Electrons and positrons annihilate each other from T = 5 GK to T = 1 
GK. Through annihilation and scattering they can transfer some energy to neutrinos, which in turn can alter 
the neutron-proton conversion reaction rates and ultimately 4 He production. Figure (3b) shows the change 
AA4 He of the mass fraction of 4 He for the same ranges of r\ and as in Figure (3a) . AA4 He imitates A"4 He , 
but not exactly. The change also remains within an order of 10 -4 . The author will determine if the change 
is equivalent to a constant shift of Figure (3a) to the left, and if the coding for the neutrino heating effect is 
correctly represented in the code. 
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